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We apply Schrodinger-functional techniques to the SU(2) lattice gauge theory with Nf = 2 flavors 
of fermions in the adjoint representation. Our use of hypercubic smearing enables us to work at 
stronger couplings than did previous studies, before encountering a critical point and a bulk phase 
boundary. Measurement of the running coupling constant gives evidence of an infrared fixed point 
where !/<?« = 0.20(4)(3). At the fixed point, we find a mass anomalous dimension 7m(<?*) = 0.31(6). 



PACS numbers: 11.15.Ha, ll.10.Hi, 12.60.Nz 



I. INTRODUCTION 

Recently, there has been an upsurge in interest in ap- 
plying the techniques of lattice gauge theory to theories 
with different numbers of colors, different numbers of fla- 
vors, and fermionic representations other than the fun- 
damental [H-Q. The immediate phenomenological appli- 
cation of these studies is to candidate theories of physics 
beyond the Standard Model, in which new gauge dynam- 
ics allow the replacement of the fundamental Higgs boson 
by a bound state of new fermionic degrees of freedom — 
so-called "technicolor" theories |6|. More generally, it is 
an interesting question in quantum field theory: How do 
systems of gauge fields coupled to fermions behave? 

Many studies attempt to characterize the behavior of 
a theory by computing a suitably defined running cou- 
pling constant. Perturbation theory 0, Q gives us a first 
approach and a menu of possibilities to be confronted by 
numerical studies. The two-loop beta function is 
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where, for an SU(iV) gauge theory with Nf flavors of 
fermions in representation R, 



bi = ^C 2 (G)-^N f T(R), 
34 

b 2 = -[C 2 {G)] 2 ~N f T{R) 
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Here C2(R) is the value of the quadratic Casimir oper- 
ator in representation R [where G denotes the adjoint 
representation, so C 2 {G) — N], while T(R) is the con- 
ventional trace normalization. Three possibilities for the 
behavior of the massless theory are: 

• triviality, b\ < in our conventions. 

• QCD-like physics, f3(g 2 ) < for all g 2 , meaning 
bi , 62 > in Eq. (TJJ . This is presumed to be associ- 
ated with confinement and chiral symmetry break- 
ing, as in ordinary QCD. 



• a fixed point at g = g*, such that f3(g 2 ) = 0. For 
b\ > this comes about if 62 < and then g* is an 
infrared-attractive fixed point (IRFP). Correlation 
functions decay algebraically at large distance, and 
there is no confinement, no chiral symmetry break- 
ing, and indeed no particle spectrum. For a given 
representation R, the domain of (N c ,Nf) where an 
IRFP exists is called the "conformal window." 

Technicolor demands a theory with QCD-like physics. 
Walking technicolor requires a QCD-like theory just out- 
side the conformal window. 

In a massive theory, the running coupling <?(/i) is aug- 
mented by the running fcrmion mass m(fi). In addition 
to the beta function, one considers the anomalous di- 
mension 7 m of the mass operator ipip. It determines the 
running of the mass parameter according to 



dm(fi) / o\ / \ 

M — j = -7m(9 MW- 



(1) In lowest order in perturbation theory, 
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In the massless theories used for technicolor, j m gov- 
erns the running of the condensate Oipip). It is thus an 
important diagnostic for realistic "extended" technicolor 
models, especially those based on walking. Phenomeno- 
logical constraints on such models require it to have a 
large, nonperturbative value. If the massless theory is 
conformal, on the other hand, then near 
correlation length £ scales as 



m q = the 



(6) 



where y m — 1 + "f m (g*) is the leading relevant exponent 
of the system [j| [l(| • (Here m q is the mass defined by 
the axial Ward identity — see below.) 

Perturbation theory cannot make reliable predictions 
for properties of systems that evolve to become strongly 
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interacting at long distances. To address questions such 
as "where is the conformal window" and "what is the 
spectrum" (for a confining theory) or "what are the crit- 
ical exponents" (for a conformal theory) requires nonper- 
turbative methods. For us, this methodology is simula- 
tion of the lattice-regulated theory. 

We present here a study of the SU(2) gauge the- 
ory coupled to two flavors of Dirac fermions in the ad- 
joint representation of the gauge group. Following the 
introduction of this theory as a technicolor candidate 
[Til H^ |. several groups have performed numerical sim- 
ulation s fl~3l- 23 . Studies of the lattice theory's phase di- 
agram [TaBll and spectrum [HM [U, QJ, indicate 
that its weak-coupling phase is quite different from that 
of SU(2) or SU(3) gauge theories coupled to a small num- 
ber of fundamental-representation fermions. Applica- 
tions of the Schrodinger functional (SF) method [13, 20 1 
and Monte Carlo renormalization group (MCRG) 23] 



have indicated that this theory has an IRFP. The relia- 
bility of this prediction is limited, however, by the lattice 
discretization used. The phase diagram obtained shows 
a first-order transition in strong coupling which termi- 
nates rather close to the claimed location of the IRFP. 
One consequence is that the value of g*, if it exists, is 
poorly determined. 

The main goal of our work was to determine whether 
this theory has an infrared fixed point. We have calcu- 
lated the beta function for the SF coupling g 2 by carrying 
out simulations on a number of different lattice volumes. 
We are able to present strong, even definite evidence for 
an IRFP. The strength of our calculation comes from 
adoption of an improved fermion action that incorpo- 
rates normalized hypercubic smearing, "nHYP fat links" 
[24l [25j . The fat-link action has been shown to effect 
a dramatic reduction of lattice artifacts when used for 
QCD simulations. As we found for the SU(3) theory with 
sextet fermions [26], this action moves the end point of 
the first order transition farther into strong coupling and 
allows us to examine a much larger range of the SF cou- 
pling without encountering it. The work of Hietanen et 
al. 1171 was limited to the range 1/g 2 > 0.22 while the 
estimate for g 2 was in the range 2.0 to 3.2 or 1/g 2 = 0.3 
to 0.5. We reach 1/g 2 ~ 0.07 to 0.10, depending on the 
volume, at the strongest bare coupling. (Here the beta 
function has changed sign, so the largest running cou- 
pling is obtained on the smallest volume.) We observe 
directly an IRFP at 1/g 2 = 0.20(4) (3) where the first er- 
ror is statistical and the second is systematic. This is a 
slightly weaker coupling than the two-loop perturbativc 
value of g 2 — 7.9 or 1/g 2 = 0.126. With this result, we 
determine the value of the mass anomalous dimension at 
the IRFP, 7 m (#*) = 0.31(6); here the bulk of the error 
reflects the uncertainty in g*. 

The outline of the paper is as follows: In Sec.|TT]we re- 
view our lattice action and the techniques we use to mea- 
sure the beta function and 7 TO . In Sec. Mil we present the 
phase diagram of the lattice theory. Sections IIVI and |V] 
contain our results for the running coupling constant and 



mass anomalous dimension. We discuss our results in the 
context of the literature in Sec. IVII 



II. METHOD 

The lattice calculations in our SU(2) study, following 
the Schrodinger functional method, were carried out in 
the same way as in our recent study of the SU(3) gauge 
theory with sextet fermions. We refer the reader to 
Ref. [26fl for a detailed presentation. In Sees. Ill Al and lll 51 
we give a short recapitulation to make this paper self- 
contained. In Sec. Ill CI we describe at some length the 
extraction of the beta function from fits to the running 
coupling. The method used here is special to a theory 
with a slowly running coupling. Our method for calcu- 
lating j m , again similar to that in our SU(3) work, is 
presented in Sec. Ill Dl 



A. Lattice action and simulation 

We study the SU(2) gauge theory coupled to two fla- 
vors of dynamical fermions in the adjoint representation 
of the color gauge group. The lattice action is given by 
the single-plaquette gauge action and the Wilson fermion 
action with added clover term [23]. The gauge con- 
nections in the fermion action employ the differentiablc 
hypercubic smeared link of Ref. [24|, from which the 
adjoint-representation gauge connection for the fermion 
operator is constructed. [The adaptation of the original 
SU(3) construction to SU(2) is trivial.] The smearing 
parameters for the links are the same as in Ref. [24] : 
Q'i = 0.75, a2 = 0.6, a3 = 0.3. The parameters that 
are inputs to the simulation are the bare gauge coupling 
P = 4/<7q and the fermion hopping parameter k, related 
to the bare mass mo by k = (8 + 2too) _1 . Tests of non- 
perturbative improvement [28-30] indicate that we can 
safely set the clover coefficient to its tree-level value of 
unity. 

The molecular dynamics integration is accelerated with 
an additional heavy pseudo-fermion field as suggested by 
Hasenbusch [3l| , multiple time scales [32| , and a second- 
order Omelyan integrator [33| . Lattice sizes range from 
6 4 to 16 4 sites. 

We study the massless theory by fixing k = k c {0), the 
point at which the quark mass m q vanishes for each /3. 
We define m q using the unimproved axial Ward identity 



(AWI), 



d t ]T (A Q (x, t)O a ) = 2m q (P a (x, t)O a ) . (7) 



where the axial current Af, 



^l^l5(T a /2)ip, the pseu- 
doscalar density P a = V'75(T a /2)i/), and O a could be any 
source. We follow the usual SF procedure and take the 
source to be the gauge-invariant wall source at t = a as 
in Eq. I|19p below. The correlation functions in Eq. ([7]) 
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are then measured at t = L/2, the midpoint of the lat- 
tice. The derivative is taken as the symmetric difference, 
d»f(x) = [f(x + fia) - f{x - Aa)]/(2o). 

On a finite lattice, the quark mass m q generally de- 
pends on the lattice size L as well as on (j3,n). As we 
show below, the dependence on L is quite weak except 
at the strongest couplings. We generally defined k c by 
demanding m q = on a 12 4 lattice. Since the L de- 
pendence becomes significant at strong coupling (see the 
Appendix), we also carried out a complete SF calcula- 
tion at a shifted n for j3 — 1.4. As will be seen below, 
the results turned out to be insensitive to this shift. 



B. Schrodinger functional and the running 
coupling 



TABLE I: Summary of simulation runs for obtaining the 
Schrodinger functional coupling g 2 at the bare couplings 
(/3, k c ), for the lattice sizes L used in this study. Trajecto- 
ries were of unit length. 



p 






trajectories 




L — 6a 


L = 8a 


L = 12a 


L = 16a 


3.0 


0.12682 


16.2K 


32.2K 


30K 


19.5K 


2.5 


0.1276 


16.2K 


32. 2K 


40.6K 


27.2K 


2.453 


0.12766 


16.2K 




16.2K 




2.445 


0.12769 


16.2K 


16.2K 






2.0 


0.12967 


16.2K 


32.2K 


41.6K 


27.2K 


1.985 


0.12279 


16.2K 




16.2K 




1.97 


0.12991 


16.2K 


16.2K 






1.75 


0.13216 


16.2K 


32.2K 


32. 3K 


41.8K 


1.5 


0.13617 


16.2K 


32.2K 


46.2K 


32. 3K 


1.4 


0.13824 


16.2K 


32.2K 


43K 


33.2K 



The Schrodinger functional (SF) [34- 38] is an imple- 
mentation of the background field method that is espe- 
cially suited for lattice calculations. It involves perform- 
ing simulations in a finite volume of linear dimension L, 
while imposing fixed boundary conditions on the gauge 
field. The classical field that minimizes the Yang-Mills 
action subject to these boundary conditions is a back- 
ground color-electric field. By construction the only dis- 
tance scale that characterizes the background field is L, 
so the n-loop effective action T = — log Z gives the run- 
ning coupling via 



T = g(L)- 2 S- 



YMi 



(8) 



where SyM i s t he classical action of the background field. 
When r is calculated non-perturbatively, Eq. ([8]) gives a 
non-perturbative definition of the running coupling at 
scale L. In a simulation, the coupling constant is deter- 
mined through the differentiation of Eq. ((8} with respect 
to some parameter r\ in the boundary conditions. This is 
an observable quantity, 



dr) 
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(9) 



The quantity K is just a number [35|. With bound- 
ary fields as described in Ref. 37], it takes the value 
K = — 12-7T in the infinite- volume limit; we use this value 
for finite volume as well, since the corrections are nu- 
merically small. We also impose twisted spatial bound- 
ary conditions on the fermion fields, following Ref. [36j j . 
4>(x + L) = exp(i9)ip(x), with 9 — tt/5 on all three axes 

The observable Q is quite noisy and requires long sim- 
ulation runs, as shown in Table HI 



C. Extraction of the beta function 

By calculating the inverse running coupling on lattices 
of size L and sL, we obtain the discrete beta function 
(DBF) 



B(u,s) 



1 



1 



1 



gi(sL) g 2 {L)' 
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(10) 



The usual beta function refers to infinitesimal scale 
changes. We define the beta function for the inverse cou- 
pling as 



/3(Vff 2 ) 



d(i/g 2 ) 

dlogL 



2/3( 3 2 )/ff 4 = 2u 2 /3(l/ U ). (11) 



Hence, 



logs 



,L dV 
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u+B(u,s) du , 



P[u>) 



(12) 



While the literature is careful to distinguish between 
the DBF and the usual beta function, we remark that 
in our theory the DBF's we can measure are, to high 
accuracy, just proportional to the beta function itself. 
This occurs for two reasons. First, our coupling runs 
slowly: We are, after all, near the (anticipated) bottom 
of the conformal window. Second, due to the cost of 
simulations, monitoring the volume dependence of the 
running coupling at fixed bare parameters is practical 
only for s < 2. 

If the beta function changes little in the course of in- 
tegrating Eq. (fT^j) . then the reseated DBF, defined as 



R(u, s) 



B(u,s) 
logs 



(13) 



will be approximately equal to the beta function /3(u). 
At one-loop order R( lS> (u,s) = — 2&i/167r 2 , a constant 
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[compare Eqs. (QJ and ([lip] . The situation at the next 
order is illustrated in Fig. [T] The figure shows the two- 
loop result, 



R (2) (u,s) = - 



2bi 

16tt 2 " 
log [1 



16tt 2 6i 
(26i/16tt 2 ) 



_1 log si 



logs 



(14) 



for the rescaled DBF for scale factor s = 2, 4, 8, com- 
pared to the one-loop and two- loop beta functions. 1 The 
rescaled DBF for s = 2 is hardly distinguishable from the 
beta function. 

There are two lessons to be drawn from Fig. [1] If 
the actual DBF resembles the two-loop result, we can 
combine the rescaled DBF's for many scale factors s onto 
a single plot to give a good approximation to the actual 
beta function. Furthermore, since any value of s < 2 is 
as good as another, we can combine the couplings for all 
lattice volumes studied to extract the beta function via 
a fit. 

We do not have to rely on perturbation theory when 
going beyond the approximation of a constant rescaled 
DBF. Expand the beta function f3(u) linearly near some 
fiducial value m, 



du 
dlogs 



j3{u 1 ) + {u-u 1 )p'{u 1 ) = B + B 1 u. (15) 



Upon integration, 

I \ at ^exp(Bilogs) - 1 
u(s) -in =p{ux) . (16) 

When the product Bi log s is small the exponential can be 
expanded, and we again observe that [u(s) — u\]/ log s = 
R(u, s) gives the beta function. 

This discussion suggests that we can fit the running 
coupling from all volumes at fixed bare coupling to 



l/g 2 (sL ) = c + a logs + c 2 (logs) 2 + 



(17) 



where L$ is a fixed reference volume. We treat all the 
parameters cq, c\, c%, ■ ■ ■ , as independent, having in mind 
that Eq. (|15p is in itself just an approximation. If the 
terms nonlinear in log s are small, the slope C\ gives the 
reduced DBF directly, according to Eq. (flU)) , We do this 
below. The success of this analysis shows again that the 
beta function in this theory is small in the region studied; 
we use other fits to estimate the systematic error. 

The discussion so far has ignored discretization errors. 
The usual analysis found in the literature presents the 
beta function only after an extrapolation to (a/L) — > 0. 
This is done by collecting data for DBF's on different 
volumes, and at the same value of g 2 . We can do this by 




ce o - 



1 In this model, b\ = 2 and 62 



-40. 



u = l/g 



FIG. 1: Rescaled discrete beta function, calculated in two 
loops for various scaling factors s. Also shown are the one- 
and two-loop beta functions; the rescaled DBF approaches 
the two- loop beta function when s — >• 1. Top to bottom, the 
curves are in the order shown in the legend. 



performing simulations at one value of the bare coupling 
f3 on our largest lattices (16 4 volumes), then moving to 
smaller volumes, matching l/g 2 (L) at fiducial volumes 
and measuring \/g 2 (sL) on appropriate larger volumes. 
For example, with s = 2 we can look at L = 16 and 8 
at f3 ~ 2.5 (for example) compared to L = 12 and 6 at 
slightly offset j3 values. We must move along the k c line 
as we do this. We have attempted this at two values of 
the SF coupling. We found (see below) that the shifts in 
(3 were not well determined with the statistics available. 
This is, in fact, a disadvantage of having such a small 
beta function. 

Such an analysis, however, is needlessly complicated 
for the question we set out to answer, namely, is there 
an IRFP, and, if so, where is it? We can test various 
hypotheses for the dependence of the running coupling 
on the lattice volume at any fixed bare coupling. Equa- 
tion (TIT)) assumes that this dependence reflects contin- 
uum physics only; in addition, we will try fit functions 
that test if the volume dependence can be explained by 
the anticipated form of discretization errors, i.e., pow- 
ers of a/L. We can also test for the presence of lattice 
artifacts by varying the data sets kept in the fit. If the 
results agree on the existence and location of a zero of 
the DBF, then we can claim to have found an IRFP. 



D. Mass anomalous dimension 

The volume dependence of the renormalization fac- 
tor Zp of the isovector pseudoscalar density P a = 
i ) lb{ T ' 1 /2)4> gives the mass anomalous dimension 7 m . 
(The pseudoscalar density is related by a chiral rotation 
to which is the object of interest.) It is computed 
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from two correlators via [2CJ, I3QM4 lj 



(18) 



" /p(£/2)' 

/p is the propagator from the t — boundary to a point 
pseudoscalar operator at time xq, 



fp(x ) 



It, 



I — T° 

d 3 yd 3 z (ip(x )-f 5 Y^( x o) 



x C(y)75yC(^) 



(19) 



We take xq = L/2. In the expression, £ and £ are gauge- 
invariant wall sources at t — a, i. e., one lattice layer away 
from the t = boundary. The /i factor is the boundary- 
to-boundary correlator, which cancels the normalization 
of the wall source. Explicitly, it is 



fi 



—y 

3L 6 



d 3 ud 3 vd 3 yd 3 z 



C(«b6yC» 



x C(2/)75yC(z) 



(20) 



where £' and £' are wall sources at i = L — a. 

We use the same boundary conditions for the calcula- 
tion of Zp as for the simulations that generate the data 
for the SF coupling. This makes its computation para- 
sitic on the SF runs. 

The constant c allows imposing a volume-independent 
normalization condition in the weak-coupling limit. Since 
we will only need ratios of values of Zp to find j m , the 
overall normalization is irrelevant. We set c = l/v2 in 
tabulating Zp below. 

We extract the anomalous dimension of tptp from the 
change in Zp [Eq. (|18|)] between systems rescaled as L 
sL. The (continuum) mass step scaling function [20|, |3S 
S3 is 



a P (v,s) 



Zp(sL) 



Z P (L) 



9 2 W 



(21) 



It is related to the mass anomalous dimension via 



dp{v, s) = exp 



j lm (g 2 (tL)) 



(22) 



Because the SF coupling g 2 (L) runs so slowly, Eq. 
is well approximated by 



(23) 



We can therefore combine many sL values collected at the 
same bare parameter values into one fit function giving 

7m, 



log Zp (L) = — 7 m log L + const. 



(24) 



This fitting procedure parallels keeping only the cq and 
ci terms in Eq. ((T7|) . As in the case of the DBF, we 
can look for subleading continuum corrections and/or for 
lattice artifacts by modifying the fit functions or the data 
set kept in the fit. 



0.15 



0.14 



0.13 



O 1st order transition, 6 

4 

E — H 1st order transition, 8 

G — © K , 12 4 

c 

A 2nd order transition, 6 4 
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FIG. 2: Phase diagram in the /3-K plane, determined with 
Schrodinger functional boundary conditions. The k c line ends 
at a point (/3*,/t*) presumed to be a critical point, where 
the indicated first-order boundary ends as well. The triangle 
marks a point on the second-order phase boundary between 
confined and deconfined phases. This phase boundary is pre- 
sumed to continue all the way to k = 0. Its upper endpoint 
may coincide with (/?*, k*). 



III. PHASE DIAGRAM 

Our determination of the phase diagram in the (/3, k) 
plane is preliminary. The diagram is qualitatively con- 
sistent with that given in Refs. [ljl [l6| . but of course it 
is quantitatively different because of the fat-link action. 

In order to measure the SF coupling, we had to map 
out the k c line. As discussed above, we determine k c by 
demanding m q — at fixed /3; this is possible only for 
sufficiently weak coupling (large /3 — see Fig. ^ . The k c 
line meets at (3 = f3* a line of first order transition at 
which the AWI quark mass jumps discontinuously from 
a positive to a negative value; this makes it impossible 
to define n c for /?</?*. On the first-order line, the 
discontinuity in m q , like that in the plaquette, varies with 
j3 and appears to vanish at fi* (see Fig. [3]); this makes the 
meeting of the two lines a critical point. This is similar 
to what was reported in Refs. [la, Il6l|. 

In finite volume, another line of transitions separat- 
ing the strong-coupling confining phase from a decon- 
fined phase begins at or near the meeting point and 
runs out towards k — 0. If the spatial volume were 
to be made large, this would be the finite-temperature 
confinement transition. The adjoint fermions leave the 
global Z(2) center symmetry unbroken, so the finite- 
temperature transition can be an Ising-like second order 



6 




m 



0.15 



0.4 


' ' ' i 1 '\m ' 




0.3 




— 


0.2 






- 
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-0.2 
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B-B 12 4 




-0.4 


i 





0.13 



0.14 
K 



0.15 



FIG. 3: Fixed-/3 scans of the average plaquette (left) and the AWI mass (right) on three different volumes. From right to left, 
the successive groups are for /3 = 1.2, 1.25, 1.3, 1.35, 1.4, and 1.5. For /3 < 1.35 there is a transition that strengthens with 
increased volume. 



transition like that of the pure gauge theory. We have 
investigated this transition at one value of k below the 
critical point and found it to be a continuous transition 
for volume 6 4 (Fig.HJ, quite different from the jumps seen 
in the fixed-/? scans of Fig. [3J The simplest scenario is 
to suppose a second-order phase boundary that stretches 
from k = to the critical point at (/3*, k*). 



IV. RUNNING GAUGE COUPLING 

Our SF calculations were performed along the k c line. 
A summary of the data collected is shown in Table Q] 
The measured SF couplings are tabulated in Table [TT] 
and plotted (for some values of (3) in Fig. [5] The log- 
arithmic variation of 1/g with L is characteristic of a 
slowly-running coupling. The transition from positive to 
negative slope as the bare coupling (3 decreases is our first 
piece of evidence for the existence of an IRFP. 



0.12 



0.1 - 



<\L\> o.08 



0.06 



0.04 



i 1 i 1 i 1 r 



_L 



n i E 



J i I i L 



_L 



1.1 1.2 1.3 1.4 1.5 1.6 
P 



A. The DBF for s = 2 

The data can be combined in various ways. The most 
direct is to plot the s = 2 DBF for each of two values of 
L/a, 6 and 8. This is shown in Fig. HJ it is essentially 
a comparison of the DBF between two lattice spacings. 
The data for L/a = 6 and 8 give a picture of the DBF 
with only weak dependence on the lattice spacing. Re- 
markably, the numerical result tracks the two-loop DBF 
rather closely. We will estimate g*, the location of the 



FIG. 4: Scan of the Polyakov loop average (\L\) at k — 0.135 
on volume 6 4 . 



zero, in our more extensive analysis below. 

At /3 = 1.4 we found that k c was strongly dependent on 
volume. We therefore calculated the SF coupling at the 
value of k c appropriate to L — 12a and then at a shifted 
value of k. We present a comparison of the two cases 
in the Appendix. The shifted k value gives a separate 
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TABLE II: Schrodinger functional couplings l/g 2 from this study. 



p 








Vff 2 




L — 6a 


L = 8a 


L = 12a 


L = 16a 


3.0 


0.12682 


0.5846(27) 


0.5771(26) 


0.5708(37) 


0.5690(66) 


2.5 


0.1276 


0.4417(27) 


0.4378(25) 


0.4273(34) 


0.4268(50) 


2.453 


12766 


0.4284(24) 




4178f50) 




2.445 


0.12769 


0.4305(27) 


0.4170(37) 
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FIG. 5: SF coupling l/g 2 vs. L/a (plotted on a logarithmic 
scale) for the four volumes studied, at (from the top) /3 = 3.0, 
2.5, 2.0, 1.75, 1.5, and 1.4. The lines through the data points 
are fits to the data at each /3 of the form l/g 2 (L) = a + 
61og(L/a). The dotted line has the slope 2b\ / (16n 2 ) as given 
by the lowest-order beta function, Eq. ([3]). 



«3 




0.4 

u = l/g 2 (6 4 or 8 4 ) 



FIG. 6: Discrete beta function, Eq. (|10[) . for scale factor s = 2 
as a function of l/g 2 measured on the smaller lattice. Two 
values of L are shown for the smaller lattice, L = 6 and 8. 
The pairs of points correspond to calculations carried out at 
(right to left) (3 = 3.0, 2.5, 2.0, 1.75, 1.5, and 1.4, as shown in 
Fig. plus points calculated at a shifted k value at /3 = 1.4. 
The dashed line is the two- loop result. The shifted k value is 
discussed in the Appendix. 



determination of the DBF at (3 = 1.4, which is included 
in Fig. |51 One sees that the shift in n leads to a change 
in the DBF that is less than the statistical error bar. 



B. Finding g* and estimating systematic error 

As can be seen from Fig. [SJ the spread of values of l/g 2 
for different volumes at fixed bare coupling j3 is typically 
much smaller than the difference between values obtained 
at different /3's. This motivates us to analyze the data in 



two stages. In the first stage, the data set at each /? is 
treated as an independent fitting problem. The outcome 
is the value of the DBF at some reference value of l/g 2 
reachable at that /3. In the second stage, the DBF's from 
all f3's are combined to obtain an estimate of the 
location of the IRFP. The variety of fits studied gives us 
a handle on the systematic error in g*. 

The data plotted in Fig. [5] are evidently linear in 
log (L/a) for fixed /3. To study deviations from linear- 
ity we try fitting four different functions to u = 1/ g 2 (L): 
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u = a + blogx, 
u = a + b log x + c(loga;) 2 , 
u = a + b log x + c/.x, 
u = a + b log x + c/x 2 . 



(25a) 
(25b) 
(25c) 
(25d) 



In these formulas, x = L/8a. The simplest fit, Eq. (|25al) . 
assumes linearity in log(L/et) and no discretization errors. 
The results of this fit at each bare coupling were plotted 
in Fig. [5] In Eq. (I25b[) we have added the next term 
from Eq. (fT7[) , reflecting subleading continuum running. 
In the last two fit functions we include instead a term 
that accounts for discretization errors. In Eq. (|25cl) we 
assume that the leading lattice artifacts are linear in a/L, 
whereas in Eq. (|25b[) they are assumed to be quadratic. 

In the continuum limit, a/L 0, fits (|25c[) and (|25d|) 
both reduce to the simplest fit, Eq. (|25a|l . The fit pa- 
rameter a is thus interpreted as l/g 2 (L = 8a), while b is 
the estimate R(g 2 ) for the beta function /3(l/g 2 ) at the 
same L. The parameters a and b have a similar inter- 
pretation for fit (|25b|) as well. In particular, since f3 for 
l/g 2 (L — 8a) is du/d\ogs evaluated at s = 1, this is 
again b. 

We plot the fit parameter b, which gives the beta func- 
tion, versus a, the inverse coupling, for each fit type in 
Fig. [7J The squares show results of fits using all four 
volumes, L/a = 6, 8, 12, 16. These fits have one degree 
of freedom, except for the simplest fit (|25al) which has 
two degrees of freedom. We have included in panel (a) 
another fit to Eq. (I25ap where only the three largest vol- 
umes (L/a — 8, 12, 16) are kept. This fit has one degree 
of freedom. 

For almost all beta values, all of these fits produce 
X 2 /dof near unity; the exception is j3 = 1.5, where all 
fits give x 2 — 6- The five fit types give results that are 
consistent with each other as well as with the two-loop 
beta function. We remark that all fits have good x 2 /dof 
at the strongest bare coupling, /? = 1.4, leaving no doubt 
that the beta function has crossed zero. 

Our next task is the determination of g*, the value of 
the running coupling where the DBF vanishes. For each 
fit type, we can locate the zero by a linear fit to the points 
in its figure. In all cases, we get good % 2 after dropping 
the points at the smallest ((3 = 1.4) and largest (/? = 3.0) 
couplings. The final results are shown in Table IIIII One 
observes that all the results are mutually consistent. 

The differences among the estimates quoted in Ta- 
ble |TlT] for 1/g 2 reflect our systematic uncertainties, to 
which we now turn. Unlike in QCD, where all sources of 
systematic error are well under control, in the case of a 
(nearly) conformal theory the systematic error is poorly 
understood. Indeed, an important conclusion from Fig. [7] 
is that our data do not allow us to sort out discretization 
errors from subleading continuum corrections. With this 
in mind, we estimate the systematic error by keeping a 
subset of the five fit types that represents both options. 



First, we obviously keep the simplest, linear fit (|25a[) 
on all four volumes. To account for the possibility of 
continuum corrections, we include fit (|25bp . 

It remains to select a fit that represents the possible 
discretization errors. Here we face a difficulty. As can be 
seen in Fig. [7J and in Table [Till fits (I23c|) and (|25d|) . which 
both include a term for lattice artifacts, give significantly 
larger error bars than the other fits. This shows that our 
data do not resolve log a: from 1/x 2 (and, even more so, 
from 1/x). We stress that our clover action with the 
nHYP links generally shows much smaller discretization 
errors than the simple Wilson action. Nonetheless, since 
we keep the clover coefficient at its tree-level value of 1, 
some residual linear dependence on a/L could survive. 

Luckily, we have yet another fit type that is sensitive 
to the discretization errors and, at the same time, pro- 
duces much tighter uncertainties in g* . This is the simple 
fit (|25a|) in which we drop the smallest volume, L = 6a. 
Since the smallest volume necessarily contains the largest 
discretization errors, dropping it must give us a result 
that is closer to the continuum limit. The advantage 
of dropping the smallest volume in the linear fit over 
fits (|25cl) and (|25dl) is that there is no need to postulate 
anything about the concrete form of the discretization 
errors; in particular, we do not have to assume anything 
about the relative size of a/L and (a/L) 2 errors. 

Disregarding fits (|25c[) and (|25dl) while keeping the 
other three, we finally conclude that 



0.20(4)(3), 



(26) 



where the first error is statistical and the second is sys- 
tematic, representing the spread of the mean values of 
the three selected fit types. It follows that 



5 n+ 2 - 7 

°- u -1.3 i 



(27) 



where we have combined the systematic and statistical 
errors linearly. 2 

The derivative of the beta function at the fixed point 
is a universal quantity. In our linear fits, reported in 
Table [TTT1 this is just the fit parameter B . We conclude 

^ = -0.08(3). (28) 

This time the error is entirely statistical. The three fit 
types we have kept produce essentially the same result, 
so the systematic error is negligible. This translates into 
an exponent 



Vg 



df3(g 2 ) 



dg 2 



= -~B =0.040(15). (29) 



2 Examination of our graphs shows that most of the systematic 
error is due to one data point for 1/g 2 , namely that at /3 = 2.0 
for L = 16a. This is responsible for the high point in the 8 — > 16 
DBF plotted in Fig. [6] at u ~ 0.3. 
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FIG. 7: Estimate R(g 2 ) for the beta function /3(l/p 2 ) as extracted from the fits listed in Eq. I)25[). Each panel shows results 
of a different fit, as described more fully in the text: (a) Squares show fits to a + blogx using all four volumes; circles, using 
L/a = 8, 12, 16 only, (b) Fits to a + blogx + c(logx) 2 . (c) Fits to a + blogx + c/x. (d) Fits to a + blogx + c/x 2 . Plotted 
curves are the one-loop (dashed line) and two-loop (dotted line) beta functions. 

TABLE III: Linear fits to the beta functions plotted in Fig. [T) Each beta function (resulting from the fitting procedure listed 
in the table) is a fit to R(g 2 ) — A + B/g 2 , and the fit gives an estimate of the zero of R at g*. 



fit type giving beta function 



.4 



B 



i/g 2 



(a) a + blogx, all volumes 

(a) a + blogx, L/a = 8, 12, 16 

(b) a + b log x + c(log x) 2 

(c) a + blogx + c/x 

(d) a + b log x + c/x 2 



0.015(5) 

0.019(7) 

0.015(8) 

0.027(34) 

0.020(18) 



-0.08(2) 

-0.08(3) 

-0.09(3) 

-0.07(14) 

-0.07(7) 



0.20(3) 

0.23(4) 

0.17(4) 

0.40(41) 

0.31(13) 
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The positive sign indicates infrared irrelevancy of the 
gauge coupling. 



C. Attempt to disentangle continuum running 
from lattice artifacts 

Here we describe an attempt at a more traditional SF 
analysis, as described in Sec. [TTJ As we have seen above, 
data taken at the same value of bare parameters, but 
at several values of L, will show a combination of true 
running and lattice artifacts. In principle, these effects 
can be separated: One adjusts the bare parameters in 
the simulation to match the SF coupling on two different- 
sized small lattices, and then increases the lattice sizes 
by the same scale factor s. A comparison of the DBF's 
obtained will show only lattice artifacts. This differs from 
our plot of the DBF for s = 2, Fig. [6j in that here we 
calculate B(u,s) for fixed u = 1/g 2 with two different 
lattice spacings; in the plot we kept (J3, k) fixed, and 
hence the lattice spacing, between lattices with L = 6a 
and L = 8a. It basically amounts to using lattice data at 
several couplings as a substitute for our fitting functions 
Ea.l25clorl25dl 

We attempted such a match near two values of the bare 
coupling, (3 = 2.0 and 2.5. In Figs. [5] and we calculate 
B(u, 2) using L — 8a and 16a at one bare parameter 
value (3, and then use L' = 6a' and 12a' at /?'; the relation 
between /3 and ft is u(L' = 6a', (3') = u(L = 8a, /3). We 
may then declare that V = L in physical units. We can 
also do the same exercise for s — 4/3 by starting with 
L = 12a and 16a at j3, matching to V = 6a' and 8a' at 
(3'. The data for 1/g 2 are shown in Fig. [8] It is apparent 
that matching the couplings, for example u(L' = 6a', P') 
and u(L — 8a, (3), can only be accomplished within large 
error bars. 

Fig.[9]shows the rescaled DBF R(u, s) for scale factors 
s = 4/3 and 2. We plot R against (a/L) 2 , and attempt 
to fit the data to a linear dependence in (a/L) 2 . We 
can do it separately for each value of s, or, following 
the discussion in Sec. HH we can fit the rescaled DBF to 
a common line. The results of those fits are shown in 
the two figures. The fits of all the data give x 2 = 0.4 
and 3.2 for two degrees of freedom at j3 — 2.0 and 2.5, 
respectively. 

At P = 2.5, or 1/g 2 ~ 0.42, the continuum- 
extrapolated rescaled DBF R(u, s) is 0.013(27) for s — 
4/3, -0.008(22) for s = 2, and 0.002(17) for the com- 
bined fit. The two loop result is —0.018, and our numer- 
ical result is consistent with it. 

At = 2.0, or 1/g 2 ~ 0.27, the continuum- 
extrapolated R(u,s) is 0.031(22) for s = 4/3, 0.027(17) 
for s = 2, and 0.024(14) for the combined fit. This is 1.7<r 
away from zero. The two-loop result is —0.014, about 2a 
away. The fits with many L's using Eqs. (|25c|) - (l25dl) . 
which include ansatze for discretization errors, also pro- 
duced a positive DBF with a large uncertainty at (3 = 2.0. 
The SF coupling is right on the edge of the value we quote 



TABLE IV: Values of Zp, the pseudoscalar renormalization 
constant 



p 






Z P 




L = 6a 


L = 8a 


L = 12a 


L = 16a 


3.0 


0.12682 


0.511(1) 


0.490(1) 


0.462(2) 


0.441(2) 


2.5 


0.1276 


0.484(1) 


0.456(1) 


0.423(1) 


0.402(1) 


2.0 


0.12967 


0.427(1) 


0.394(1) 


0.354(1) 


0.332(1) 


1.75 


0.13216 


0.373(1) 


0.338(1) 


0.301(2) 


0.277(1) 


1.5 


0.13617 


0.293(1) 


0.261(1) 


0.227(1) 


0.203(1) 


1.4 


0.13824 


0.250(3) 


0.219(2) 


0.188(1) 


0.173(2) 



in Eq. (|2"6"|) for 1/g 2 from our analysis of many couplings. 

Notice that the analysis in this subsection underesti- 
mates the error in the DBF since we have not included 
the uncertainty in j3' that arises from matching the SF 
couplings. To make a definitive determination of g 2 us- 
ing the extrapolation method would require repeating 
it at many couplings with significantly better statistics. 
We conclude that this method is no better than fits to 
Eqs. (|25cp - (|25d|) . We have already argued that the other 
fits, Eqs. P5a |l -([25b ll . give more reliable results. 



V. MASS ANOMALOUS DIMENSION 

Finally we turn to the mass anomalous dimension j m . 
Our analysis will parallel that of Sec. IIVBI with the dif- 
ference that log Zp(L) replaces l/g 2 (L). 

Our data for Zp(L) are shown in Table HVl and plotted 
in Fig. QUI The lines represent fits of the form of Eq. (|24|) 
to the data at each bare parameter value. The nearly 
straight-line behavior of the data is indicative of a slowly 
running coupling. 

Fig. QT] shows the values of 7 TO extracted using 
Eqs. (l2~Tj) and (|2"3")l from pairs of lattices with scale factor 
s — 2. We have plotted the data as a function of the SF 
coupling g 2 measured on the smaller volume. There does 
not seem to be a great deal of difference between results 
from the two pairs, meaning that there is not much of a 
shift with lattice spacing. The rightmost points are from 
the strongest coupling, j3 = 1.4. As is clear in our various 
determinations of the DBF, the SF coupling does run sig- 
nificantly at this value of (3. Note again that the L = 6a 
coupling is larger than the L = 8a coupling, indicative of 
the positive DBF. 

The main feature of Fig. [IT] is that the measured 
anomalous dimension j m at first follows closely the one- 
loop curve, but beyond g 2 ~ 4 it flattens out. 

As we did for the running coupling, we can fit 
log Zp(L) to various functional forms. Again we be- 
gin with the simple linear behavior of Eq. (|2"4")l ; we pro- 
ceed to add corrections with the aim of testing whether 
deviations reflect continuum corrections or lattice arti- 
facts. We use the same fit functions, given in Eq. (|25l) . 
with the change that u in these equations now stands 
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FIG. 8: 1/g 2 near /3 = 2.0 (left) and /3 = 2.5 (right). Points at = 2.0 and /3 = 2.5 are slightly displaced for clarity. 




FIG. 9: Discrete beta functions near j3 = 2.0 (left) and j3 = 2.5 (right), where the smaller lattices' g 2 values are matched by 
slightly varying the bare parameters. The points near (a/L) 2 = are the results of fits, described in the text: the circle is the 
extrapolation of R(u, 4/3), the square is the extrapolation of R(u, 2), and the diamond is a fit to all four points. Some points 
have been displaced slightly for clarity. 



for log Zp(L). The reasoning that led to the identifica- 
tion of the fit parameter b with R{g 2 ) in Sec. IIVBI now 
leads to its identification with — 7m(g 2 )- The results are 
shown in Fig. [T2J In all cases we plot 7 m as a function 
of g 2 (L = 8a), where the latter was obtained from the 



parallel fit type in Sec. IIVBI 

Our data for Z p are far more precise than the data for 
1/g 2 given above. This is reflected in the much smaller 
error bars in Fig. compared to Fig. [7] At the same 
time, this leads to rather high values of x 2 /dof (in the 
range 3-4) in the linear fits (I%4"1) . whether we use three 
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FIG. 10: Pseudoscalar renormalization constant Zp. From 
the top, data are from f3 = 3.0, 2.5, 2.0, 1.75, 1.5 and 1.4. 
Lines are fits to log Zp(L) = — 7 m log(L/a) + const for each 

1 . i i i i i i i i i i i i i i i i i i i i i 



or four volumes. The \ 2 is good for the other three fits, 
again excepting j3 = 1.5 which gives x /dof ~ 6. 

Finally we estimate 7m (ff*), the mass anomalous di- 
mension at the IRFP. The uncertainty in 7 m (<?*) must 
reflect the uncertainties in the determination of g* itself. 
In Fig. [13] we have plotted again our results for the linear 
fit to all four volumes. The horizontal bar at the top of 
the figure indicates our result for g^, Eq. (|2"7|). We see 
that the right and left ends of the error bar almost co- 
incide with the points corresponding to j3 = 2.0 (on the 
weak-coupling side) and (3 = 1.5 (on the strong-coupling 
side). This observation renders unnecessary any interpo- 
lation of the curve. In view of the monotonicity of 7 m (g) , 
we simply take the values of j m at these two couplings 
to mark off the uncertainty of 7m (<?*), concluding that 

lm{9*) = 0.31(6). (30) 

Roughly 5/6 of the error comes from the uncertainty in 
<7*. The small statistical and systematic errors of 7 m 
itself are responsible for the rest. Based on the same 
reasoning as in Sec. IIVB| we took the spread of values 
obtained using the two linear fits [panel (a) of Fig. [TJ] 
and fit (|25b[) [panel (b)] as a measure of the systematic 
uncertainty in 7 TO . 

Fig. [T3] also shows a comparison of our data to those 
of Bursa et al. [2(| , where we have applied to their data 
the linear fit (|24p. (Their published graphs include an 
additional large systematic uncertainty.) Our data lie 
slightly below theirs. 
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FIG. 11: Anomalous dimension 7 m calculated from ap(v,s — 
2) using Eq. (|23p from pairs of lattices: L/a = 6 — >• 12 
(squares) and L/a = 8 — > 16 (diamonds). It is plotted as 
a function of the SF coupling g 2 on the smaller lattice size. 
Bare couplings range from f3 = 3 on the left to 1.4 on the 
right, as in Fig. 1101 The line is the lowest-order perturbative 
result. 



VI. DISCUSSION 

Our simulations show that SU(2) gauge theory coupled 
to two flavors of adjoint fermions lies inside the confor- 
mal window. Using Schrodinger functional techniques, 
we have determined its IRFP and measured two expo- 
nents, the mass anomalous dimension 7 m and the ex- 
ponent y g of the (irrelevant) gauge coupling. Both are 
small. 

Because the massless theory is conformal, not confin- 
ing, it is not a candidate for a technicolor theory. Even 
if the theory were to be deformed in a way that destroys 
the IRFP, the small anomalous dimension would spell 
trouble for its application as a technicolor theory. 

We are aware of several estimates of and 7 m (g*) in 
the literature. The two-loop zero of the beta function 
occurs at gl — 7.9 which is right on the edge of our 
quoted range. Our central value is at weaker coupling, 
Vl = 5.0. 

The location of the IRFP is scheme dependent, of 
course. The previous lattice study with the best data, 
which also used the SF coupling to find a fixed point, 
is that of Hietanen et al. [17| ■ Their published analysis 
combines data from all Us and bare parameter values 
into a single fitting function. They quote fitted values 
of gl in the range 2.0-3.2, or l/gl in the range 0.3-0.5. 
We are in mild disagreement with their results. Their 
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FIG. 12: 7 m from fits of logZ p plotted as a function of g 2 (L — 8a): (a) Squares show fits to a — j m logs, where x = L/8a, 
using all four volumes; circles, using L/a — 8, 12, 16 only, (b) Fits to a — 7,„ log a; + c(loga;) 2 . (c) Fits to a — j m log a; + c/x. 
(d) Fits to a — 7 m log x + c/x 2 . The line in each case is the one-loop result. 



lattice action uses unimproved Wilson fermions with no 
clover term, and hence is susceptible to 0(a) discretiza- 
tion errors. The line of first order transitions also extends 
farther into weaker coupling in their case; they have no 
data below 1/g 2 = 0.22. 

As seen in Fig. [T31 our results for j m ( n 2 ) agree with 
the earlier determination of Bursa et al. [20J, where we 
overlap in couplings. We also agree with the MCRG re- 



sult of Catterall et al. [2^|, j m {g*) = 0.49(13), given its 
larger uncertainty. Del Debbio et al. [U, [22| use the scal- 
ing of spectral observables with fermion mass to compute 
7 m at g 2 ~ 3, finding small values [0.05-0.20 and 0.22(6) 
respectively]. These results also lie on the lowest-order 
perturbative curve. 

We find that jmig 2 ) deviates from the lowest-order 
perturbative formula for g 2 > 4. As we saw for SU(3) 
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FIG. 13: Mass anomalous dimension 7 m (s ) from the linear 
fit (|24|) . which are the squares in Fig. I12f a~). The horizontal 
bar at the top marks our result ([23 for gl. The crosses are 
the data of Bursa, et al. [2(|, analyzed with the same linear 
fit. Again the diagonal line is the perturbative result. 



with sextet fermions, the numerical results lie below the 
curve. 

Now let us consider analytic predictions for 7 m (<?*). 
Ryttov and Shrock [i^] have an extensive tabulation of 
perturbative results up to four loops in the MS scheme 
(see also [43[ )■ Recall that only anomalous dimensions 
measured at fixed points are scheme-independent. Per- 
turbative predictions of course depend on the order of 
perturbation theory. They tabulate 7m (ff*) = 0.820, 
0.543, 0.500 for two, three, and four loops. They also give 
a prediction based on solving a Schwinger-Dyson equa- 
tion, jm(g*) = 0.653. The Ryttov-Sannino all-orders 
beta function Q gives 7m(g*) = 0.75. Pica and Sannino 
offer another all-orders beta function [45[ which gives 
7m(<?*) = 0.46. All these numbers are too high to agree 
with our result. 

There is an extensive literature attempting to relate 
the location of the bottom of the conformal window to a 
large value for j m (see 0, |46[ and references therein) . 
Our results indicate that models using the SU(2) gauge 
group and adjoint fermions are not relevant to that lit- 
erature: the parameter space accessible to exploration 
is too granular. Our 7 m is small. One might want to 
decrease Nf in the hope that 7 m would grow as one ap- 
proaches the bottom of the conformal window. This is 
clearly not possible here: The Nf = 1 theory has b 2 > 
and is probably confining. 

Finally, we point out that the use of an improved ac- 
tion, which smooth the gauge field fluctuations as seen 



by the fermions, made this project feasible. It shifts 
the location of the strong-coupling transition deeper into 
strong coupling than could be achieved with the simple 
Wilson action, allowing us access to the strong coupling 
side of the IRFP. 
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Appendix: Determining n c 

For each value of the bare coupling /3 we determined 
k c by demanding m q = for L = 12a. We calculated m q 
in a series of short runs over a range of k's, eventually 
confirming m q — in a run of several hundred trajecto- 
ries. When we used this value of k c in the much longer 
runs that yielded the SF coupling g 2 , the error bar on 
m q naturally shrank and thus the final result for m q is 
always a little bit off zero. 

We used the same values of (/?, k c ) for different lattice 
volumes, so as to keep the lattice spacing fixed for the 
SF calculation. In the weak-coupling region, the volume 
dependence of k c is weak, as seen in Fig. 1141 In particu- 
lar, m q on the 16 4 lattice is small enough that we can say 
that m q ~ in the infinite- volume limit. This is not the 
case, however, at stronger couplings (Fig. [15]) . Of par- 
ticular concern is the large value of m q for the strongest 
coupling, /3 = 1.4, on the largest lattice, L = 16a; the 
value k c = 0.13284 is perfectly adequate for L — 12a but 
not for L = 16a (see Table E}. 

In principle, one could fix k c by demanding that m q — > 
in the infinite volume limit, which can be done by fitting 
to finite- volume results at various values of k. As is clear 
from Fig. [T51 the finite-lattice corrections to m q are pro- 
portional to (a/L) 2 , as may be expected for the clover 
action. This procedure would involve lengthy calcula- 
tions, however, since the error bars shown in the figure 
are only attainable with the statistics of a full SF sim- 
ulation. We decided instead to check on the sensitivity 
of our results to a shift in k of the order that might be 
required by this procedure. 

Table fVl and Fig. [TBI show the values of m q a calculated 
at the nominal k c — 0.13824 and at a shifted value k = 
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TABLE V: AWI masses m q a at /3 = 1.4 for the two k values considered. 



p 


K 




m q a 








L = 6a 


L = 8a L — 12a 


L = 16a 


1.4 


0.13824(= k c ) 


0.092(2) 


0.035(1) 0.0013(4) 


-0.0092(3) 


1.4 


0.1381 


0.109(3) 


0.053(2) 0.0128(6) 


0.0014(3) 




TABLE VI: Schrodinj 


;er functional coupling 


5s l/g 2 at /3 = 1.4 for the two k values. 




/3 


K 












L — 6a 


L = 8a L = 12a 


L = 16a 


1.4 


0.13824(= k c ) 


0.0655(29) 


0.0790(21) 0.0950(27) 


0.1035(34) 


1.4 


0.1381 


0.0610(21) 


0.0722(24) 0.0891(34) 


0.0926(34) 



0.006 



0.004 



0.002 



m a 

q 



-0.002 



-0.004 



-0.006 



T 



G-OP=3 
B-H P=2.5 

P = 2.0 
H — h (3=1.75 




_l I L 

16 12 8 
Lla 



0.1381 that gives a small value of m q a for the largest 
lattice L = 16a. Table PVTl lists the SF couplings l/g 2 for 
all four volumes at both values of k. While the shift in k 
does induce a significant and systematic shift in l/g 2 , the 
change in the DBF is much smaller than the statistical 
error, as is seen in Fig.[Bl A similar result obtains for Zp 
(Table IVlTf : the ratios ap are unaffected and thus the 
estimate of 7 m is unaffected as well. We conclude that 
our results are insensitive to such shifts in n, even at the 
strongest coupling studied. Note that we have not made 
use of any of the data at ft = 1.4 in determining or 
lm{g*)- 



FIG. 14: Volume dependence of the AWI mass m q at fixed p, 
where k was fixed by a preliminary determination of m q — 
on the 12 4 lattice. 
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